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Abstract 

The Coulomb-gauge vector potential of a uniformly moving point charge is obtained by 
calculating the gauge function for the transformation between the Lorenz and Coulomb gauges. 
The expression obtained for the difference between the vector potentials in the two gauges is 
shown to satisfy a Poisson equation to which the inhomogeneous wave equation for this quantity 
can be reduced. The right-hand side of the Poisson equation involves an important but easily 
overlooked delta-function term that arises from a second-order partial derivative of the Coulomb 
potential of a point charge. 



1. Introduction 

Gauge invariance is an important property of electrodynamics. Notwithstanding the thor- 
ough attention it has received in textbooks [1-4], the topic of gauge invariance seems to 
be still very much alive. Recently, articles have appeared on the resolution of apparent 
causality problems in the Coulomb gauge [5], the transformation from the Lorenz gauge to 
the Coulomb and some other gauges [6] , the Coulomb-gauge vector potential in terms of the 
magnetic field [7] , and the historical development of the whole concept of gauge invariance 
[8]. A discordant voice in this is a claim by Onoochin [9] that the electric field of a uni- 
formly moving point charge comes out differently when it is calculated in the Lorenz and 
Coulomb gauges, which that author takes as evidence that the two gauges are not physically 
equivalent. 

Curiously, it is difficult to find in the literature an explicit expression for the Coulomb- 
gauge vector potential of a uniformly moving charge. To the present author's knowledge, 
such an expression has appeared only in a short comment [10], where it is obtained directly 
from the well-known electric field and Coulomb-gauge scalar potential of the charge. The 
formulae of Jackson [6] for the transformation between the Lorenz and Coulomb gauges thus 
have come timely to provide an analytical check on the claim of Onoochin. The purpose 
of the present paper is to clear up the problem of the Coulomb-gauge vector potential 
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of a uniformly moving charge using Jackson's results as well as the method employed by 
Onoochin, in which the difference between the Coulomb- and Lorenz-gauge vector potentials 
in this problem is to be found by solving a Poisson equation. 

In section 2, we obtain the difference between the potentials in the Lorenz and Coulomb 
gauges for a uniformly moving point charge using the formalism of Jackson's paper [6], 
which guarantees that the two gauges will yield the same electric and magnetic fields. In 
section 3, we demonstrate that the expression obtained for the difference between the vector 
potentials in the two gauges satisfies the inhomogeneous wave equation for this difference, 
and in section 4 we show where and why the procedure used by Onoochin for solving that 
wave equation by reducing it to a Poisson equation went wrong. It will turn out that the 
right-hand side of this equation must include an easily overlooked delta-function term to 
yield the correct solution. Concluding remarks are made in section 5, and an appendix 
contains some calculational details. 

2. Transformation from the Lorenz gauge to the Coulomb gauge 

The scalar potential V and the vector potential A of a uniformly moving point charge in 
the Lorenz gauge, defined by the condition 'V'A + dA/cdt — 0, are well known. For a point 
charge q moving with a constant velocity v = vx, the Lorenz-gauge potentials Vl and 
are given by (we shall use the Gaussian system of units): 

where 7 = (1 — j c?')"'^!'^ and, for simplicity, the charge is assumed to pass through the 
origin x=y=z=^ at time t = (see, e.g., [3], section 19-3). In the Coulomb gauge, defined 
by the condition V-A = 0, the scalar potential Vc of the charge takes a particularly simple 
form, 

^c(r-, t) = ,! „ „ (2) 



\j{x — vty + y'^ + 



as the scalar potential in this gaiigc is exactly the same as that of the instantaneous Coulomb 
interaction of electrostatics. On the other hand, the Coulomb-gauge vector potential Ac is 
the retarded solution to a relatively complicated inhomogeneous wave equation 

4:71 dVr 

□Ac = -qvS(r - vt) + (3) 

c cot 

where □ = — d"^ /c^dt^ is the d'Alembertian operator and qv6{r — vt) is the point-charge 
current density of the present problem (see, e.g., [1], section 6.3). 

The gauge invariance of electrodynamics implies that there is a gauge function xc that 
connects the Coulomb- and Lorenz-gauge potentials by 

Vc^Vi^-^ Ac = Al + Vxc (4) 

which ensures that the Lorenz-gauge and Coulomb-gauge potentials will yield the same 
electric and magnetic fields. This is because the fields are generated from the potentials via 
the prescription 
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E{r,t) = -VV{r,t) - B{r,t) = VxA{r,t) (5) 

and thus any electric-field and magnetic-field differences that could arise from the use of the 
different gauges are guaranteed to vanish: 

Vx(Ac-Al) = V X Vxc = 0. (7) 

Nevertheless, it should be instructive to demonstrate explicitly that it is indeed the case also 
in the present problem by finding the requisite gauge function. Before we turn to this task, 
we give here for completeness the fields that the prescription (5) yields with the Lorenz-gauge 
potentials Vi^ and Al of equation (1): 

^(^>') = ^ l,.(,-ll[r4'+.y^ B{r,t) = lvxE(r.t). (8) 

The same expressions for the electric and magnetic fields can be obtained also by Lorentz- 
transforming the electrostatic Coulomb field of the charge from its rest frame to the 'labo- 
ratory' frame. 

Jackson [6] has derived the following integral expression for the gauge function xc in 
terms of the charge density: 

Xc{r, t)^-cj dV'-^ ^""^^ drp{r', t - r) (9) 

where R = \r — r'\ (a gauge function is defined to within an arbitrary additive constant, 
which we omit here). For a point charge q moving with a constant velocity v along the 
X-axis, the charge density is 

p(r, t) = q5{r - vt) = q6{x - vt)6{y)6{z). (10) 

This gives 

f-R/c pR/c 



I ' dTp{r',t-T) = qd{y')d{z') f " dr S[x' - v{t - t)] 
Jo Jo 

= ^^5{y')5{z') j^''' dr 5[t - {t - x'/v)] 
= ^^S{y')S{z'){e[R/c -{t- x'/v)] - e[-{t - x'/v)]} 
= ^S{y')6{z')[Q{x' - xo) - Q{x' - vt)] (11) 
where Q{x) is the Heaviside step function and 



X — vt-\ — \ (x — vt)'^ -\ — -{y"^ + z^) 
c V 7 



^/l-v^/t 



The gauge function (9) with the charge density (10) is thus 
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c [°° , , Qix' - xq) - Qix' - vt) 
Xc{r, t) = -q- / dx' \ °; ^ ^ 



c 

= -q- 

V 



vt 



C 

Q- 

V 



^0 ^J{x- x'Y + y'^ + z'^ 

(13) 



. ^ X — vt . ^ X — Xq 

arcsmn — arcsmn 

VZ/ + V ?/ + 



Let us first check that the gauge function (13) yields the established difference Vc — Vl 
between the scalar potentials in the two gauges. Using the identity 



l-J{x-xoy + y'' + z^ ^-{x-vt) + J {x - vty + \{y^ + z') (14) 
7^ * c y 7 

which follows from (12), to simplify the result of the partial differentiation —dxc/cdt, we 
obtain 

V^-V^^JJ^^ 'I - 'I (15) 

cdt ^^x-vty + y^ + z^ ^{x-vty + ^{y^ + z^) 

This is indeed the correct result [see equations (1) and (2)]. 

Calculating the x-component of the difference Ac — Al between the vector potentials 
in the two gauges is now very simple because dxc/dx — —{l/v)dxc/dt on account of the 
dependence of the gauge function (13) on the variables x and t only through the combination 
X — vt. We thus have 

dxc cdxc c 
Ac. -^^^ = ^---^ = -(^c - V^). (16) 

The y- and z-components of the difference Ac — Ai^ — Vxc a-re obtained by performing 
direct differentiations in a similar way to that of calculating the value (15) for the difference 
Vc-Vl, yielding 

^o, - ^L, = ^ = -'-'4i4(yc - vo (17) 

ay V y^ + z^ 



Ao.-A,, = ^ = -'-'4^(Va-V,). (18) 

oz V y^ + z^ 



These components have no singularities (they vanish at x — vt — y — z — 0). As Ai,y — 
Alz — 0, equations (17) and (18) also give the Coulomb-gauge components Acy and Acz 
themselves, respectively. 

It is instructive to perform the Lorentz transformation of the Coulomb-gauge four- 
potential (Vc, Ac) from the 'laboratory' frame, where the charge moves with the constant 
velocity v — vx, to its rest frame (the primes denote the rest- frame quantities): 

V^ = ^{Vc-vAcJc) (19) 
A'c,^^{Ac.-vVc/c) A'cy = Acy A'c, = Acz. (20) 
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We note first that with the x-component Ac x of the vector potential in the Coulomb gauge 
given by 

Acx = -{Vc-Vjj') (21) 

V 

which follows from the expression (16) for the difference Aqx ^ ^Lrr and the fact that 
^Lx — {'v/c)Vi„ equation (19) gives the rest-frame scalar potential as 

V' = —V = = (22) 

where the Lorentz transformation of the coordinates, x' = 7(x — vt), y' = y, z' = z, is 
performed on the right-hand side. The rest-frame scalar potential Vq is simply that of a 
point charge q in electrostatics. We note also that because the Coulomb-gauge condition 
V'A = is not Lorentz invariant, the rest-frame vector potential A'q is not divergenceless 
in the rest-frame variables x', y', z'; the potentials Vq, Aq are therefore no longer those of 
the Coulomb gauge. However, a direct calculation shows that the vector potential A'q is 
irrotational, V'xAq = 0, which expresses the fact that there is no magnetic field in the 
rest frame. Moreover, the vector potential A'q is independent of the rest-frame time t', and 
thus the electric field in the rest frame is given only by E' — — V'Vq, which yields correctly 
the electrostatic Coulomb field of a charge at rest. 

3. Inhomogeneous wave equation for the vector-potential difference 
The difference Aq — Ai^ must satisfy the inhomogeneous wave equation 

□ (Ac - Al) = (23) 
which is obtained by subtracting the wave equation for Al, 

4:71 

□Al = qv5(r - vt) (24) 

c 

(see [1], section 6.3) from the wave equation (3) for Aq- It is straightforward to show that 
the x-component Acx — Ai^^ = {c/v){Vc — Vl) of the difference Ac — Al indeed satisfies 
the x-component of the inhomogeneous wave equation (23): 



□ 



V 



cdtdx 

The fact that equation (25) holds true follows directly from the wave equations 



(25) 



and 

□Vl = -47rqS{x - vt)5{y)S{z). (27) 
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The wave equation (26) in turn holds true because of the facts that the d'Alembertian 
□ = - 9Vc^5^^ Wc = -Anq6{x-vt)5{y)5{z), and dVc/dt = -vdVc/dx; the wave 
equation (27) embodies the fact that the Lorenz-gauge scalar potential is the retarded 
solution of the inhomogeneous wave equation with the right-hand side —4:7rq5{x—vt)S{y)S{z). 

It is also straightforward to show that the y-component (17) and z-component (18) of 
the difference Ac — Ai^ satisfy the inhomogeneous wave equation (23) by using the identity 

oifg) = gDf + fUg + 2V/.V5 - 2^^ (28) 

and equation (25) in the evaluation of the requisite derivatives; the singularities of the 
functions y{x—vt)/{y^+z^) and z{x—vt)/{y'^+z'^) y — z — Q cannot introduce any delta- 
function terms in the d'Alembertians of the components (17) and (18) as the latter are 
functions with no singularities. 

4. Poisson's equation for the vector-potential difference 

Onoochin [9] attempts to solve the x-component of the inhomogeneous wave equation (23) 
for the difference Aq — Al directly by reducing it to a Poisson equation, which is a method 
based on the fact that the space and time partial derivatives are not independent when the 
source term is moving uniformly (see [3], section 19-3). In the present case, the dependence 
on the variables x and t is only through the combination x — vt, and thus the x-component 
of the difference Ac — Ai^ can be written as Acx — Ai^x — f{x~vt,y,z), where the function 
f{x — vt, y, z) satisfies an inhomogeneous wave equation 



nf{x - vt,y,z) 



d'Vc 



cdtdx 

c [[X — vty + + z^Yi^ 3c 

that can be cast as a Poisson equation on the substitutions d"^ jdx^ — d'^ /(?dt^ — 9^/7^ 
and X — vt — x/7, where 7 = (1 — v^/c^)~^/^: 

' +#7 + 1^ /(x/7,y,^)--g- /.V7, L. +Q^lS{x)S{y)S{z). (30) 



\dx^ dy'^ dz'^ ) ' ' c (x^/7^ + Z/^ + -2^)^^^ 3c 

The delta-function term on the right-hand side of equation (29) arises from the fact that 
d^Vc/cdtdx = —{v/c)d'^Vc/dx^ and the delta-function identity [11] 

-^d^jd{r) (31) 



dxidxj r 3 

where r = {xj + x^ + Xg)^/^, Sij is the Kronecker delta symbol and S{r) = S{xi)S{x2)S{xs) 
is the three-dimensional delta function. Strictly speaking, the first term on the right-hand 
side of equation (31) should be understood as lima^o{SxiXj — r'^5ij)/{r'^ + a^)^/^, where the 
limit is to be performed after an Yi? integration with a well-behaved 'test' function. This 
limit is automatically implemented when the integration is done in spherical coordinates 
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and the integration over the angular variables is performed first (without any additional 
transformation of variables that would shift the origin, of course) [11,12]. 

Unlike the well-known delta- function identity V^(l/r) = — 47r5(r), the identity (31) 
is needed only relatively rarely in electromagnetism, an example being the calculation of 
the fields of electric and magnetic dipoles from their potentials (see [11]; and [1], equations 
(4.20) and (5.64)). The identity V^(l/r) = — 47r5(r) obviously requires that the second-order 
partial derivatives d'^r^^/dxf, z = 1, 2, 3 are given by expressions like those of equation (31); 
the precise form of these can be seen most easily to follow from the limit a — > of 

d 1 3xiXj T Sij Oj Sij (32) 



as here the limit a — > of the second term on the right-hand side is a representation of 
The standard integral expression for the solution to Poisson's equation (30) is given by 



ni} poo roc roc 

47rC J-oo J-oo J—oo 



-T^Kx)5{y)Kz) 



X , ^ (33) 

^lix-x'Y + {y-y'Y + {z-z'Y 

where one recovers the original variables by putting x = 7(^ ~ '^'t)- Here, the delta-function 
term is readily integrated to yield 

12 ^12 



f=^r dv' r d^' r d.- ^^"ir-y^-z''^ ^ 

^ Anc J-oo ^ J-oo ^ J-oo (x'77^ + y" + z'^fi^ ^lix-xT + {y-y'Y + {z-z'y 



I (34) 

3c 

A direct evaluation in closed form of the three-dimensional integral in (34) does not seem 
possible. However, this integral can be evaluated for the special case y — z — and the 
result is (see Appendix): 



^ r dv' r dv' r dz' '^^^h^-y'^-^" ^ 

A-KC J-oo ^ J-oo J-oo ix"h' + y" + ^(^-^>)2 + ya + ,>2 



qvf qv 
3c\x\ 3c|a; — vt\ 



(35) 



For y — z — 0, the value of the second term on the right-hand side of (34) is exactly equal 
and opposite to the value of (35), and thus the solution (34) of the Poisson equation (30) 
vanishes at y — z — 0, as required by the result (16) that we obtained for the difference 

Onoochin's calculation omits the delta-function term that arises from the second-order 
partial derivative of the Coulomb-gauge scalar potential Vq, and he takes erroneously just 
the first term of (34) as the solution of the Poisson equation (30). According to equation (35), 
this incorrect solution is non-zero at y = 2; = 0, and, moreover, its partial time derivative 
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does not vanish there. As d{Vc — Vi,)/dx = at ?/ = z = 0, the result (35), if it were the true 
difference Aq^ — ^Lx at y = 2; = 0, would lead there to a non-zero difference between the 
Coulomb- and Lorenz-gauge x-components of the electric field. This is Onoochin's evidence 
against the equivalence of the Lorenz and Coulomb gauges. (In his paper [9] , Onoochin does 
not evaluate the integral in equation (34), but he has communicated to the present author 
a calculation of its value at y—z—0 that differs from the value given by equation (35) by a 
factor of 3 [13].) 

It is not difficult to show directly that 

f--( ■ - . (36) 



is the vanishing- at-infinity solution of the Poisson equation (30). Using the relation (31) 
and the differentiation rule d'^g{ax,y, z)/dx^ = a'^d'^g{u,y, z)/du^\u=ax, we have 



^ 1 _ 1 

dx^ y'xVT' + y^ + z^ " 



- ^7^(x)^(2/)^(^) 



(37) 



1 V - xVt' - 47r ,,,,,,,,, 

92 1 _ 2z2-xV7'-y' 47r 



^^2/^2+^2 + ^2 (xVt^ + y'^ + Z^yi"^ 3 

which gives 

/ 92 92 92 \ 



^lKx)Ky)^{z) (39) 



\dx^ 9y2 dz^J ^^2/^2 + ^2 + 

= - (1 - ^) (^;^;;v. - f + ^) ^^MHvn^)- (40) 



Using the relation V^(l/r) = — 47r5(r), we have also 



-47r75(x)(5(t/)(5(z). (41) 



Subtracting equation (41) from equation (40), multiplying the result by qc/v and using the 
fact that 1 — — v'^ /c^, we obtain 

/ 92 92 52 \ c / q 07 



9y2 a^2y y \./xV72 + ?/ + ^2 Vx^ + 1/2 + ^2 



which shows that the function / of equation (36) indeed satisfies the Poisson equation (30). 
Transforming back to the original variables through x = 7(2^ — vt), the function / becomes 
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V 



X 



^-{Vc-Vj^) (43) 



which is the value (16) for the difference Acx — M^x- 

The function / of equation (36) has to equal the integral solution (34), yielding for the 
integral in (34) a closed-form expression 



loo ^ loo ^ i-oo " (x'Vt^ + y'-" + z'-f^^ ^{x-x'Y + {y-y'f + {z-z'f 

1 [-, \ 7 



y2 



^^2/^2 + ^2 + ^2 \^ 3cy Vx' + y'^ + z 



(44) 



Direct numerical three-dimensional quadrature of this integral, with the first factor of the 
integrand regularized as mentioned above in connection with equation (31), resulted in values 
that were close to those of the closed- form expression on the right-hand side of (44) . 

5. Concluding remarks 

We found an explicit expression for the gauge function of the transformation between the 
Lorenz and Coulomb gauges for a uniformly moving point charge. The Coulomb-gauge 
potentials obtained using the gauge function arc guaranteed to yield the same electric and 
magnetic fields as the well-known Lorenz-gauge potentials of the charge. The expression 
obtained for the difference between the vector potentials in the two gauges satisfies the 
Poisson equation to which the inhomogeneous wave equation for this difference reduces after 
a transformation of the variables. However, the right-hand side of the Poisson equation 
involves a delta-function term which has to be included in the integral expression for its 
solution. 

Although gauge invariance is a foregone conclusion in a gauge-invariant theory, we believe 
that an explicit demonstration of the equivalence of the Lorenz and Coulomb gauges in the 
basic problem of a uniformly moving charge was instructive. 
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Appendix 

After the transformation x'/l = ^'^ the integral (35) can be written as {qv'j/c)I{'j, X), 
where I{'y,X) is the integral (we drop the primes on the integration variables): 

\ roo roo roo Ix"^ y^ z"^ 1 

""^^T. L loo i-oo {x^ + y^ + z^r^^^2^x-xy + y^ + z^ ^''^ 

where X = xll ^i^*^ 7 > 1 arc real parameters. 

Transforming from the Cartesian coordinates x, y, z to the spherical ones, r, ^, 0, we have 
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r^dr / sm^d^ — = 



7(7, X) = - I dr sii 
2 Jo Jo 



^72(X-rcos^)2 + r2(l -cos2 ^) 



1 

2 Jo 



dr dC ^ 



-1 



(46) 



We integrate with respect to ^ first (this is needed to implement the limit a — > of the first 
term on the right-hand side of equation (32), see [11,12]): 



se-i 



37a; A+E + 1+27 )[uj r -37 X In — — — ^ — 

^uj\X—r\-\-uj'^r—^'^X 



where 



a; = ^72-1 A= {ujh -^^'^X)\X + r\ 
The function F(r, 7, X) has the properties 



B 



(a;V + 37^X)|X-r|. 



F(r, 7, X) = F(r, 7, -^) hm F(r, 7, X) = 0. 

r— >0+ 



(47) 



(48) 



(49) 



It peaks at r = where its derivative with respect to r is discontinuous. The argument of 
the logarithm in (47) reduces to (7|X| — a;r)/(7|X|-|-a;r) for r < and to (7— a;)/ (7+0;) = 
(7 — ujY for r > \X\. 

The integration with respect to r is performed in two parts: 



hh,X) = - I drF(r,7,X) 

1 /■^'^ 3^ 
h{^,X)^-j drF{r,^,X)^ 



2 + 7^)(l + 27^) , 1 + 27^ 

+ „ .,..1 ln(7 - uj) 



2uj^\X\ 



1 + 272 



where again a; = (7^ — 1)^/^. The whole integral 1(7, X) is thus 

1 



ln(7 — uj) 



I{^,X) = h{^,X) + h{^,X) 



(50) 
(51) 

(52) 



and the integral (35) has the value {qv^ /c)I{'^,X) — qv/{?>c\X\) — qv^ / {?>c\x\) ■ 

The results of this Appendix were obtained using the software system Mathematica [14], 
and were checked by performing numerical integrations. 
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